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Abstract 
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threshold logarithms to next-to-leading order. From this (triple-differential) distribution, any other distribution 
is obtained by integration. As an application of our method, we derive simple analytical expressions for a 
few distributions, resummed to leading approximation. It is shown that the shape function can be directly 
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energy spectrum, which has a "Sudakov shoulder" , and we show how the distribution in the singular region is 
related to the shape function. We also present an improved formula for the photon spectrum in B — ► X s 7, 
which includes soft-gluon resummation and non-leading operators in the effective hamiltonian. We explicitly 
show that the same non-perturbative function — namely the shape function — controls the non-perturbative 
effects in all the distributions in the semileptonic and in the rare decay. 
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1 Introduction 



In this note we discuss factorization and threshold resummation to next-to-leading order (NLO) in semi-inclusive 
B decays: 

B -> X q + (non QCD partons), (1) 

where X q is any hadronic final state containing the light quark q = u, s coming from b fragmentation. The 
non-QCD partons are a lepton-neutrino pair or a photon. 

In sec. 2 we consider the most general distribution in the decay 

B —> X u + I + v. (2) 

A general expression for the triple differential distribution is presented, from which any other distribution is 
obtained by integration. The process (||) is characterized by three energy scales: 

m B , Q, m x , (3) 

where^ 

Q = 2E X . (4) 
Factorization is achieved with two different steps, according to the method developed in [Q: 



1. We factorize heavy mass effects by taking the limit 

tjiq — ► oo, Q, mx — ► const. 



(5) 



The distributions contain logarithms of the heavy mass, logms/Q, related to the non-conservation of the 
transition currents in the static limit for the beauty quark This step is discussed in sec. 2.1. 



2. We factorize the large infrared logarithms that appear in the semi-inclusive region, 

a£V k {z) (0<fc<2n-l), 

by taking the limit 



where 



m 



x 



We have defined: 



V k (z) 



Q 2 



log" (1 - z) 
1 - z 



(6) 



(7) 



(8) 



(9) 



J + 



Plus-distributions are defined as usual as P (z), = P (z) — 5 (1 — z) dx P (x) . This step is discussed in 
sec. 2.2. 



1 The factor 2 multiplying Ex is inserted for convenience. 
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The basic idea of our approach is to use the kinematical variables 



w = — (0 < to < 2) and z, (10) 
m B 



the latter replacing the generally used one 



u=l-^f. (11) 

m B 



In other words, we consider (two times) the hadronic energy Q as the hard scale of the process, rather than the 
b quark mass rrif, (or, equivalently, the B-meson mass). The reason for this choice is that the infrared structure 
of the decay is not modified by the limit ([5J): this implies that tub cannot be the relevant hard scale; the 
heavy mass acts only as an energy resevoir for the light partons in the final state and has not any fundamental 
dynamical meaning. Using the variable z instead of u largely simplifies the logarithmic structure. In the latter 
variable, the decay distributions contain, at one loop, terms of the form ||: 

X>i(u), \ogwV {u), V (u), log 2 w<S(l-u), logwS(l-u), (12) 

while in the variable z one has only terms of the form: 

V 1 {z), V (z), log k;<5(1-z). (13) 

In the latter case, the two different logarithmic structures basically decouple. 

We explicitly show that a unique function / (z) , depending only on z, factorizes the long-distance effects in 
any distribution in the decay (^). The perturbative expansion of / (z) reads: 

f(z) =5{l-z)-a s A 1 T>\ (z) + a s B x V (z) + 0(a 2 s ), (14) 

where the constants multiplying the distributions are given by: 

7T 4 7T 

and Cf = {N 2 — l)/(2iV c ) = 4/3. This function is related, by a short-distance coefficient function C, to the 
shape function in the effective theory j|, ||, defined as 

tp (k + ) = (B (v) \ht S (k+ - iD + ) h v \B («)), (16) 

where 

9 

K ^ --f. (17) 

C is defined by the relation ^, |^] : 

tp (k+; Q) QCD = Jc(k + - k' + ; Q, //) <p (k' + ;n) dk' + , (18) 

where we have defined 

V {k + -Q) QCD = ~f(z) (19) 

and (i < Q is the ultraviolet cut-off or renormalization point of the effective theory. The shape function factorizes 
the long-distance effects — both perturbative and non-perturbative — in the process up to the scale \x. The 
relation between the QCD and the effective theory variable is 
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The relevance of the variable z to describe the long-distance effects — perturbative and non-perturbative — in 
(0) is also supported by the following argument. From dimensional analysis, the shape function is of the form 

V (k+;fx) = -^(^±). (21) 



A* V A 4 

The coefficient function in eq. (|18|) does not contains large logarithms of the ratio Q/ ' fj, if we choose 

M ~0(Q), (22) 

implying that 

h 



" D2 



(23) 



The last equation implies that long-distance effects are described by a function of the variable z, Q.E.D.. Let 
us stress that f(z) "includes" the shape function but does not coincide with the latter, because it also contains 
some short-distance effects (soft gluons with energies bewteen fi and Q and hard collinear contributions) || [jj . 

The shape function was originally derived with an Operator Product Expansion (OPE) and was claimed to 
be the universal non-perturbative component in (|l|) since the early papers on the subject (|, ||: we present here 
an explicit proof of this. On the perturbative side, our analysis can be considered as the resummed analogue 
of the distributions computed to O (as) in Q; some of the decay distributions presented in |jj were computed 
before in § g [fij fl. 

In the second part of sec. 2 we present simple analytical expressions for a few distributions in (|^), resummed 
to leading logarithmic accuracy. An interesting distribution from the theoretical side is the hadron energy 
spectrum. This distribution has a singularity in fixed-order perturbation theory inside the physical region, in 
the point |j| 

E X = (24) 

This phenomenon originates from the mismatch of the allowed phase space in lowest order in as :[] Ex < 
tub I 2, and in higher orders: Ex <mB- The all-order resummation eliminates the singularity and produces a 



characteristic behaviour known as the "Sudakov shoulder" |L2] . Actually, a singularity very close to point ( 24 ) 
remains even after soft-gluon resummation. The latter has a different origin: it is related to the Landau pole 
in the running coupling and is factorized with the introduction of the shape function. 

In sec. 3 we present an improved formula for the photon spectrum in the rare decay 

B -> X s + 7 (25) 

in the threshold region. We also explicitly show that the same function, the effective form factor / (z) , factorizes 
all the long-distance effects in processes (^ and (p5|). 

The distributions in @ and (^H|) can be roughly divided into two classes. The first set contains distributions 
not involving integration over the hadronic energy. A first example is the z distribution in (Q): 

d[l-z) T>i[z)-- T> (z)-\ , (26) 



Tq dz 7T 4 7T 

where To = G\ m\\V u b\ 2 / (1927T 3 ) is the total semileptonic width in Born approximation. A second example is 
the photon spectrum in the rare decay: 

1 dT RD asC F 7a s CV , , . . 

= 6(1-Xy) 2?i(x 7 )-- V (x 1 )^ , (27) 



T RD dx 1 7T 4 7T 



2 Analogous computations were performed a long time ago in the context of one-loop QED corrections to decay jllj]. 

3 Kinematics involve the decay of a massive particle into three massless particles and it is analogous to the well-known case 



where T RD is the total b — * 57 width and 

2,E 

x 7 = — 1 (0 < x~ < 1) . (28) 

771B 

In the decay (^5|), the hadronic energy is never integrated because kinematics fixes Q « raj (see sec. 3 for a 
proof). As expected, the leading terms, T>\ (z) and 2?i (x 7 ) , have the same coefficient. It is instead non trivial 
that the subleading terms also, 2?o (z) and T)q (x 7 ) , have the same coefficient, —7/4, and that these distributions 
have a perturbative expansion similar to the one for / (z) (cf. eq. (|l4|)). In general, by measuring distributions in 
this class, one can directly determine the effective form factor / (z) , or equivalently, the shape function ip (fc+) . 

The second class contains distributions in which the hadronic energy is integrated over. As examples, 
consider the hadron-mass distribution Q in (||) 



Tq du 7r 12 

or the electron spectrum || in (||), 



^ (f u) D\ (u) —— 2?o (m) + • • • , (29) 



= 1 S lo S C 1 - ^e) - 77: log (1 " X e ) + • • • , (30) 



2r dx e . 2vr ' 12 



where 



(0<a? e <l). (31) 



= — 

m B 

The distributions (|2^) and (^0|) have the same leading behaviour as in (|2^) and (p7|); the coefficient of the 
subleading terms, —31/12, is instead different because the integration over the hadronic energy affects them 
Jill . The distributions in this class are not directly related to the shape function and the extraction of the latter 
from the experimental data requires in general a deconvolution. 

Finally, in sec. 4 we present our conclusions and an outlook of future developments. 



2 Semileptonic decay 

Let us consider the hadronic tensor containing all QCD dynamics: 

W MV =J2(B\J+\X U )(X U \J^\B) {2nf5\ PB -q-p x ) 1 (32) 
x u 

where q is the momentum of the lepton-neutrino pair and (x) = u (x) 7 M (1 — 7 5 ) b (x) is the b — > u current of 
the Standard Model. The latter involves five independent form factors. We find it convenient to use a modified 
parametrization with respect to the one proposed in || : 

W^ v (p B ;Px) = 2i> 1 p x ^ n ^ Vv + n " V ^ ~ 9nvV -n- ie f _ lva 0n a v 13 ) W\ (g, w) - g^ v W 2 (<r, w) + (33) 

+v^v u W 3 (?, w) + {n ll v v + n v Vy) W 4 (?, w) + n^Ws (?, w)] , 

where 

it 

«*» = -£S- = (l; 0,0,0) (34) 
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is the velocity of the beauty meson, which we take at rest, while 

(1;0,0,-V?) 



Px 



v -p x 

is the normalized momentum of the jet, which we have taken in the minus direction. We have defined: 



(35) 



x 



(o < C < i) 



(36) 



x 



Note that 1 — q = 4(1 — z). We have inserted for convenience a factor l/(2v ■ px) multiplying all the form 
factors. Note that n ■ v = 1 and that n 2 = 1 — <r <C 1 for c < 1, i.e. n is close to the light cone in the threshold 
region. 

The form factors can be decomposed as: 
Wi{z; w; a s ) = 6 a 5(1 - z )6(l-w)+ S tl ^^- s (z) + ^^6(1 - z) Vl (w) + n ((; w) + 0(a 2 ). (37) 

7T 7T 7T 

The function s (z) contains the plus distributions, i.e. the long-distance contributions: 

7. 



s (z) = -V x (z) - -I> (z) . 



(38) 



The functions u, (w) come from virtual effects — they are proportional to 6(1 — z) — while the functions fj (£; w) 
originate from real emission only. The explicit expressions for these functions can be extracted from || and are 
reported in the following. The "virtual" functions read: 



V\ (w) = 

v 2 (w) = 

v 3 (w) = 

V4 (w) = 

v 5 (w) = 



3 w\0gW 5 7T 2 

--log^-L^l-^)-^-^-^^ 

0; 
0; 

1 



w log w 
2(1 -w) \ l-w 



1 - 2w 



2(1 -w) \l-w 



log w — 1 



(39) 



where Li 2 (z) = YlnLi z ™ l n2 (\ z \ — -0 ^ s the standard dilogarithm. The "real" functions are: 



r x (C; w) 
ri (s,w) 
r3 (?, w) 

r 4 w) 

T5 {S,w) 



w 2 (8 - w) (2 - w) 
~4 4( 
w (8 -w) 32-8w + w 2 H (?) 



w (2 - w) (8 - w) (2 - w) 



8? 



H {<;) + - \H(<;)+\og(l-z)] 

1 — z 



8 8c 16 

w(8-3w) 32 + 22w-3w 2 3w(12-w 



w 2 C + 2w(4-w) + 



32 - 8w + w 2 



g(0 

16 



id (I 5w (4 — w) 



4c 8c 2 

64 + 56u> - 7w 2 3w (12 - w) 



w 2 u>(32-5w) 3uj(12-k;) wH (?) 
1 8c + 8c 2 16 - 



8? 8c 2 
w(8 + w) 3w(12-w) ia0"(?) 



— 3w — 



2(22-3w) 3(12-io) 



8c 



16 



2 (2 -to) 3(12-uj) 



(40) 
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We have defined: 



"M^T^" <41) 

In the semi-inclusive region: 

ff(0 = -Iog(l-2)+O(l-«) ) (42) 

implying that there is no 1/ (1 — z) singularity in n for z — ► 1. The functions rj have at most a log(l — z) 
singularity for z — * 1, which gives no logarithmic enhancement after integration over z. Another property is 
that the "real" functions do not have any singularity of the form l/( 2 for £ — * 0, but at most a singularity of 
the form 1/Q (the point ? = corresponds to the final hadronic system at rest). As already noted, the Born 
term in eq. (p7|) has the additional restriction w < 1. 



2.1 Heavy mass effects 

Let us now consider the limit of infinite mass for the beauty quark, keeping the other kinematical invariants 
fixed, i.e. the limit (|^). In terms of our variables, this is: 

w -> 0, C const. (43) 

In the limit ( ^3|) the "virtual" functions behave as: 

3, 

vi -> --logw; 

v 2 -> 0; 

«3 -> 0; 

I 

«4 - 2 5 

w 5 -» 0. (44) 

The "real" functions have limits: 

ri -» ~ + |ff(0 + r ^[^(c)+log(l-*)]; 

4 2 
8 4 , . 
r 4 -» 0; 

r 5 -» 0. (45) 

Only the function v\ diverges logarithmically in this limit: all the other functions have a finite limit or vanish. 
The "real" functions are finite because real emission diagrams do not generate logarithms of the heavy mass, 
so logu> does not appear. The logarithmic divergence of v\ is associated to the term in W\ : 



Wi = 6(1 - z) 



3a s C F 
1 logw 

2 7T 6 



+ • • • . (46) 



Equation ( [Iq ) can be understood by considering the general properties of the vector and axial current containing 
a heavy and a light field. The matrix element of these current between quark states in QCD are of the form 



G 



Ml: 



and 



(u\V^\b) = u u 



{u\An\b) = u u 



a s C F f 3 

7 log m B + 

7T \ 4 



a s C F ( 3 
1 H 7 log m B 

7T \ 4 



In 75 



2tt 



-^75 r 



(47) 



(48) 



where the dots denote terms dependent on the kinematics of the external states. Substituting the curly brackets 
in the above equations in place of the currents in eq. (|32|), we recover the logarithmic term in eq. (^) for Wi, 
together with the delta- function contribution in W4. 



The appearance of the logarithm of the heavy mass in the matrix elements ( \17\ ) and (48) can be understood 
by considering the effective theory: it is related to the non-conservation of the vector and axial currents in the 
static limit 



nib — > 00. 



(49) 



The currents are multiplicatively renormalized in the effective theory |2| and their matrix elements between 
quark states read: 



{u\V v (/i) I b) = u u 
{u\A v (11) \b) = u u 



a s C F ( 3 
1 + 7 log /^ 

7T \4 

a s C F ( 3 
1 + 7log/i 

7T \ 4 



(50) 



where fj, is the renormalization point and the dots denote /i-independent terms. As expected, the logarithmic 
dependence on the renormalization point [i of the effective currents matches the dependence on the heavy mass 
mt in the matrix elements of the full QCD currents. 



2.2 Infrared factorization 



To perform factorization of infrared logarithms, the form factors are conveniently written as: 
Wi(z; w; a s ) = S a f (z; a s ) + 5(1 - z) v t (w) + n (C; w) - 6 a 6(l - z)6 (w - 1) + 0(a 2 s ). (51) 

The function / (z) , defined previously (eq. ([l4])), contains the long-distance effects. The soft-gluon resummation 
of f(z) allows us to describe the semi-inclusive region 

1 - z < 1 (52) 

and has been performed to NLO in [l^, 0); f(z) is perturbatively computable as long as 

A 



where A is the QCD scale. In the region 



1- z > 



1 - z 



Q' 



A 

Q 



(53) 



(54) 



f(z) acquires a substantial non-perturbative component related to the well-known Fermi- motion effects 
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The moments of the effective form factor, 

f N = f dzz N - 1 f(z), (55) 
Jo 

can be written as the exponential of a series of functions: 

f N = cxp [L gi {/3 a s L) + g 2 (/3 a s L) + a s g 3 {P Q a s L) H ] , (56) 

where the leading and next-to- leading functions have the simple analytical expressions |[o| [l5| [7|: 

91 W = -^ L l[(l-2A)log(l-2A)-2(l-A)log(l-A)]; 

gAX) = - M [bg(l _ 2A) _ 2bg(l _ A)] _ [lQg2(1 _ 2A) _ 21og2(1 _ A)] + 

+Al og( i_2A) + ^log(l-A). (57) 

We have defined L = logn and n = N/N 0) with N = e"^ = 0.561459 . . . and -y E = 0.577216 . . . the Euler 
constant. The first two coefficients of the /3-function are: 

_ \\C A -1n F _ 33-2n F _ 17C% - 5C A n F - 3C F n F _ 153 - 19np 

Pa ~ 12^ ~ 12tt ' lh - 24^ ~ 24^ ' (58) 

where C A = N c — 3 and n F = 3 is the number of active quark flavours. The one- loop quantities Si and C\ 
have the values: 

5 1 = -^, <* = -7— • (59) 

7T 4 7T 



The two- loop quantity A2 is given by 16 1?J : 



A2 = %K (60) 

Z7T 



where, in the MS scheme for the coupling constant, 



with T R = 1/2. 

The original form factor / (z) , in momentum space, is obtained with an inverse Mellin transform of /jv in 

/(*)=/ =7 (62) 

where the constant c is chosen so that all the singularities of /jy lie to the left of the integration contour. The 
inverse transform is usually done numerically |l8| ]. In leading order, the transform from /jv to the cumulative 
distribution 

F(z)= [ f(z')dz' (63) 

J z 

is equivalent to the simple replacement: 

logn^-logy, (64) 
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where 



m , 



The leading form factor can then be explicitly given in momentum space and reads: 



(65) 



/0; a> s )i = — exp [h (y; a s )} , 
ay 



(66) 



where 



h(y;a s ) = - [(1 + 2(3 a s log y) log (1 + 2p a s log y) - 2 (1 + /3 «S log y) log (1 + [3 a s log y)} . (67) 

It holds that h (y = 1) = 0. In agreement with the frozen coupling case (see later) we define: exp [h (y — 0)] = 0. 
The QCD coupling is evaluated at the hard scale of the process: 



as = ct s (Q 2 ) = as (w 2 m 2 B ) . 
To some approximation, we can set w = 1 in the running coupling so that 

a s ^ a s (m|) . 

Equation (pq) is simple but non-trivial, as the presence of the singularity for 



V < Vn 



exp 



1 



2(3 Q a s 



A 



(68) 



(69) 



(70) 



shows. This singularity is related to the infrared pole in the running coupling. The frozen-coupling case is 
obtained by taking the limit /3 — * on the r.h.s of eq. (pTj) and readsfl 



dy 



exp 



usCf , 2 
log y 



2tt 



asC F logy 
it y 



■ exp 



«sCf , 2 
— log y 



Wo = 0) 



(71) 



Our task is to factorize the long-distance contributions found in the hadronic tensor, i.e. the terms containing 
distributions. The simplest factorization scheme involves a minimal subtraction, so that the form factors are 
written: 

Wi (c, w; a s ) = Vi (w; as) f(z; as) + Ri (?, w; a s ) ■ (72) 

According to the above decomposition, the hadronic tensor is: 

W^u (<r, w; a s ) = Vp V (w; a s ) f(z; a s ) + R^u (<r, w\ a s ) ■ (73) 

The tensors containing the virtual effects and the real ones arc defined in a way analogous to the hadronic 
tensor: 

V lxv (w;a s ) = 2v ~ [{ n n v v + n v v n " 9^v -n- ie tiua i 3 n a v ) Vi (w) - g^ v V 2 (w) + (74) 

+v fi v u V 3 (w) + (n M tv + Vi (w) + n^rivVs (w)] , 

4 The plus regularization of the function log (1 — z) j (1 — z) in the last member can be omitted because the exponential suppresses 
the virtual contribution for ag ^ 0. 
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and 



R uv (s,w;a s ) = 2v ^—^ [{n a v v + n u v^ - g^ v v -n- ie^ va pn a v^) R x (c, w) - g^ R 2 to) + (75) 

+t^u„ i? 3 (?, to) + (n M «„ + n v v a ) i? 4 to) + 7t M n„i? 5 (?, to)] . 



The form factors have an expansion in powers of 015 



Vi (to; a s ) = 5u + ° S - v t (w) + (—) v[ (w) + 

7T \ 7T / 



and 



asCF f as\ z 

Ri k, to; as) = -5a 6 (1 — z)6 (w — 1) -\ to) + — r, (<r, to) + 

7T V 7T / 



(76) 



(77) 



All the dependence on the heavy mass, i.e. on to, is contained into the short-distance form factors: the function 
f(z; as) depends only on the ratio of the final-state kinematical variables. 



In leading approximation, the hadronic tensor has the particularly simple form: 
Wpu {z,w;a s ) l = 
with h (y; a s ) given by eq. @ 



— - — (naVu + n v Vn - guv v-n — ieu Ua pn a v ) ^- exp [h (y; a s ) 
tub w ay 



(78) 



2.3 Triple differential distribution 

Let us now consider the most general distribution in (H), which is a triple differential distribution. One has 
basically to contract the hadronic tensor with the leptonic one. A third kinematical variable is involved, which 
we choose as the electron energy: x = x e . The expression in terms of the form factors reads: 



1 d 3 T 
12Tq dxdwdz 



[z, to, x 



a s ) = /,Pi (x, to, z) Wi (to, z; as) , 



where Pi (x,w,z) are polynomials in all the kinematical variables (independent of as): 



(79) 



Pi (x, to, z 
Pi (x, to, z 

P 3 (x, w, z 
Pi (x, to, z 
P5 (x, to, z 



= [1 + x — w] [to — x — (1 — z) w 2 ~\ ; 
= ^[l- W + (l-zM; 



~[x(w-x)-(l- z)w 2 ] ; 
x (to — x) — (1 — z) w 2 ; 

— [x (tO — 2?) — (1 — z) «J 2 ] . 



(80) 



Explicitly, one has: 
1 d 3 T 



12Tr, dxdwdz 



(z,w,x;as) = 



[1 + x - to] [to - x - (1 - z) w 2 ] Wi + — [l - to + (1 - z) to 2 ] W 2 



+ [x (to — x) — (1 — z) to 2 



^W 3 + w 4 + -w 6 

4 to 



(81) 
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where x = 1 — x. 

Let us now consider the kinematical constraints; there are various cases ||. For a given electron energy, in 
the range 

0<z<l, (82) 

the range of the hadronic energy is 

x<w<l + x, (83) 

and the range of y is: 

max 0, —)<y< a ■ 84 

For a given hadronic energy iu in the range 

< to < 2, (85) 

the range of y is: 

it; — 1 



0. 



< y, (86) 



and the range of the electron energy is: 



| (1-V?)<3?<! + (87) 
For a given hadronic mass C in the range 

< «■ < 1, (88) 

the range of the hadronic energy is: 

w < — (89) 
1 + V? 

and the range of the electron energy is again given by eq. ( j37| ) . 

In the semi-inclusive region, the triple differential distribution is naturally written as: 

i d 3 r 

12T~ dxdwdz ( z > w > x > as } = C {™, x \ a s) f (z;a s ) + D (q,w,x;a s ) . (90) 
The coefficient function and the remainder function have a power series expansion in as : 



C(w,x;as) — c (w, x) + as ^ F c ( W; x ) + f — \ c'(w,x) + --- 

7T V 7T / 

D (z, w, x; as) = d (z, w, x) + ( — I d (z, w, x) H . (91) 

The function D has at most a logarithmic singularity ~ log* (1 — z) for z — > 1 and its lowest order vanishes: 
do (z, w, x) = 0. The explicit expression of the coefficient function reads: 

C (w, x; a s ) = ^2 F i ( x , w, 1) Vi [w; a s ) = Pi (x, w, 1) H P { (x, w, 1) Vi (w) H . (92) 

i=l i=l 
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Expanding in powers of a§ on both sides, the first two terms of the coefficient function read: 



c (w,x) 
c (w, x) 



Pi (x, w, 1) = (1 + x — w) (w — x) ; 

5 

^2 P * ( X , W A) Vi (w) 



(1 + x — w) (w — x) 



3 wlogw 5 7T 2 

-_] Q gt -I Ja (l-t )-__^- i - T 



+ x (w — x) 



(93) 

(94) 

log w 
2(1- w)' 



The remainder function reads: 

5 5 

D (z, w, x; a s ) = Y2 Pi ( x > w ' z ) Rl ( w ' z ' as ) + X! ^ ( x ' w ' z ) ~ ^ ( x ' w ' ^ ( w; as ^ f ( z > as ^ ' ( 95 ) 



Since 



Pi (x, w, z) - Pi (x, w, 1) = O (1 — z) , 



(96) 



one can replace the fixed-order expansion of / (z; as) in the r.h.s. of eq. ( p5| ) and neglect the virtual contribu- 
tions, i.e. the plus regularization. The one-loop contribution to the remainder function is: 



d (z, w, x) = y~] P t (x, w, z) n (w, z) + [Pi (x, w, 1) — Pi (x, w, z)] 



i=i 



log (1 - z) + 7/4 
1 - z 



(97) 



The explicit expression for d is quite long and we do not report it here. 



Equation (J9TJ) is our main result and allows computing an arbitrary distribution in the threshold region to 
NLO logarithmic accuracy. It is the generalization to a triple differential distribution of the representation for 
the resummed shape variables for very small values of the resolution parameters, such as the thrust distribution 
for 1-T« 1 



In leading order, one needs the effective form factor / (z; as) in double-logarithmic approximation, eq. (p6[), 
and the coefficient function at the Born level given above, cq (w, x) , in eq. (^3|). Replacing the expressions for 
the coefficient function and the effective form factor, the distribution reads, in leading order: 



i d 3 r 

To dxdwdy 



= 12 (1 + x — w) (w — x) f (z; as) 

= 12 (1 + x - w) (w - x) ^ ex P [h(y;as)\ , 



(98) 



with h (y) given in eq. (67) 



In next-to-leading order (NLO), one needs / (z; as) to single-logarithmic accuracy and the one-loop functions 
c(w,x) and d(z, w, x) . Since the term containing the long-distance effects depends only on z, the integrations 
over the electron energy and the hadronic energy do not touch the infrared logarithmic structure. 



2.4 Distribution in the hadron and electron energy 



Integrating over the variable y on both sides of eq. (p8|), in the range (p4|), one obtains for the leading distribution 
in the energies w and x: 



1 d 2 T 

To dxdw 



12 (1 + x — w) (w — x) < exp 



x (w — x) 



(Aw) exp 



Aw 



(99) 
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where 



Aw = w — 1. 



(100) 



Since in the term proportional to 9 (Aw) , the hadron energy is restricted to the range 1 < w < 2, according to 
the leading logarithmic approximation, we can set w 2 = 1. Introducing the neutrino energy x v = 2E 1/ /rriB in 
place of the hadron energy, satisfying 



x e + x v + w = 2, 

the distribution can be written in the more "symmetrical" form: 



1 d 2 T 

Tq dx e dx v 



V2x v (1 — x v ) < exp 



h 



(1 - X e ) (1 ~ X„) 
(2 X e 



(101) 



1 (1 — x e — x v ) exp [h (1 — x e — x v )\ > . (102) 




Fig. 1: Plot of the z distribution in leading logarithmic approximation in the range 0.75 < z < 0.95. Solid line: 
running-coupling case; dotted line: frozen-coupling case. 

Since we are interested to the semi-inclusive region 1 — x e <C 1, we can set x e = 1 whenever possible, so that: 



1 d 2 T 



Y2x v (1 — x v ) I exp 



1 - x e 
1 - x,. 



Tq dx e dx v 

In the frozen-coupling limit, the above expression reduces to: 

i d 2 r 



— 6 (1 — x e — x v ) exp [h (1 



(103) 



Tq dx e dxv 



Ylxv (1 — x u ) i exp 



O-sCf , 2 I 1 - X r- 

■log 



2?r 



1 — Xv 



— (1 — x e — x,/) exp 



■log (1 - x e - I„) 



2tt 



(104) 



The expansion to first order in as reads: 

-pr , d t ~ I2xv (1 - Xv) \ 8 (x e + x v - 1) - Q '! C ' F log 2 f 7 — — ] + 9 (1 - x e - ar„) - log 2 (1 - x e - ac„) 
lodx e aXv I 2n \1 — Xv / 2tt 

(105) 



This formula contains the same double logarithms as the one-loop distribution computed in || . 
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2.5 Distribution in the hadronic variables z and w 



Performing the integration over the electron energy in the range fl87|), the double distribution in z and w reads, 
to leading logarithmic accuracy: 



1 d 2 T 

T Q dwdy 



2w 2 (3 - 2w) f (y) = 2w 2 (3 - 2w) — exp [h (y; a s )] 



(106) 



Note the complete factorization in the variables w and y. Measuring this distribution, it is possible to determine 
in a direct way the shape function and to check the factorized structure. 




Fig. 2: Plot of the z distribution in the running coupling case for z > 0.95, i.e. for very large z. The dotted 
line represents the real part of the distribution, as the latter acquires an imaginary part after the peak. 



2.6 Distribution in z 



Integrating the above distribution over the hadronic energy in the range (p9|), 

w<l + 0(y), (107) 
one obtains the leading distribution in the (square of the) hadron mass normalized to the hadron energy: 



T dy dy 



(108) 



Note that this distribution does not coincide, beyond the double logarithm at one loop, with the hadron mass 
distribution, as discussed in the introduction (cf. eq. (p9[)). Since / (y) is proportional to the shape function via 
a short-dista nce f actor, a measure of this distribution allows a direct determination of the shape function. The 
distribution (108) is plotted in fig. 1. The region very close to z = 1 is plotted in fig. 2. The maximum of / (z) 
occurs at 



(109) 



i.e. in the non-perturbative region described by the shape function^]. The leading distribution (108) acquires 
indeed an imaginary part for z > z max . 

5 1 wish to thank S. Catani for a discussion on this point. 
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2.7 Hadron energy spectrum 



Integrating over y the distribution (10£) according to condition ( p6| ) (the upper limit on y can be taken to be 
one for simplicity's sake), one obtains the resummed hadron energy spectrum in leading approximation: 



1 dT 



= 2w 2 (3 - 2w) {1 - 6 (w - 1) cxp [h (w - 1)]} . 



(110) 



To dw 

The above distribution acquires an imaginary part — and therefore is completely unphysical — in the region 



< Aw < exp 



1 



2(3 a s (m|) 



A 

m B 



This effect is related to the Landau pole in the coupling and implies that the region 



Ex G 



tub m B 
~2~' ~Y 



0(A) 



(111) 



(112) 



is non-perturbative. This region is described by the shape function because it is related to the integration over 
y down to Aw -C 1. 



In the frozen-coupling case, one has the simple expression: 



1 dT 
To dw 



2w 2 (3 - 2w) < 1 - 6 (w - 1) exp 



2tt 



log 2 (w - 1) 



(/J - 0) . 



(113) 



The expansion to order as of the r.h.s. is in agreement with the double logarithmic approximation of the 
one- loop distribution computed in pfl: 



1 dT 
Tq dw 



2w 2 (3 - 2w) { 6 (1 - w) + 6 (w - 1) 



1 a s C F 2 7 a s C F 

log (W — l) + — log [w — 1) 



(114) 



Note the factor 7/4 in front of the single logarithm, which is characteristic of the distributions not integrated 
over the hadron energy, as discussed in the introduction. The hadron spectrum is plotted in fig. 3. Let us 
remark that, in the frozen-coupling case, the distribution is smooth in the point w = 1, because the exponential 
in eq. (113) goes to zero together with all its derivatives for w — > +1. A plot of a small neighborhood of this 
point would indeed show that the apperent cusp in fig. 3 is actually not there. 



3 Rare decay 

In this section we present an improved expression for the photon spectrum near the endpoint in the rare decay 
( |25| ) [ p0| , [2l| , |22|. The only terms in the spectrum containing large infrared logarithms to O (as) involve two 
insertions of the operator Oi (see |2^] for the definition of the operator basis). The latter is also the dominant 
one in the inclusive rate, so this operator usually is the only one considered. It is however possible to improve 
the distribution by including the other operators, as we are going to show. 

The hadronic variables z and w introduced to describe the scmilcptonic decay are redundant in this case, 
because the photon has q 2 = 0; we have: 



1 



(2-xf 



1 — x, 



w = 2 - x w 1, (115) 
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where x = x~. The variables z and x basically coincide in the semi-inclusive region while the hadron energy is 
fixed, as anticipated in the introduction. We therefore consider the x spectrum. The improved distribution we 
propose reads: 



dT Gla, 



,m iVt^fml^mlj^im,) [Q(p h ) f (x) - p' (x)] + O (a 2 s ) , 



dx 32tt 4 



where 



(116) 



r.(P) i„„ TOb 



(117) 



The "remainder" function p' (x) is the derivative of: 



p(x) 



_ a s 0&) 



£cf 0) ( Mb )cf fo) /«(!-*) 



The latter vanishes at the endpoint: 



(118) 



p (x) — > for x — ► 1, 



(119) 



because the functions /y (1 — x) all vanish in the endpoint: 



fij (1 — x) — * for x — > 1. 



(120) 



The function (a;) has at most a logarithmic singularity ~ log (1 — x) for cc — » 1. The quantities and C- 
are the leading and next-to-leading contributions to the effective coefficient functions: 



ri are complex constants depending on the ratio m c /m&. The definition of all the symbols can be found in 



(121) 




Fig. 3: Plot of the energy spectrum in leading logarithmic approximation. For w > 1, the dashed line corresponds 
to frozen coupling, the continuous line to running coupling. In the latter case, the distribution is not well defined 
close to the point w = 1 because of Landau-pole effects. 
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To avoid large logarithms in the matrix elements, one has to take: 



/i 5 = O (mb) 



(122) 



The long-distance effects are contained in the function / (x) , which receives all the non-perturbative contribu- 
tions (in leading twist) related to Fermi motion. The main point is that / (x) is the same function / (z) , which 
was introduced to factorize the long-distance effects in the semileptonic decay (cf. eq. (po|)). Therefore we 
explicitly verify the universality of the long-distance contributions in (jl]), once they are factorized by means of 
a function of z. The function p' (x) is instead a short-distance contribution, specific of the process (p5l). 



The proof of eq. (116) is trivial. The integrated distribution in the upper part of the photon spectrum, 

T up (x) = j —dx', (123) 
given in eq. (30) of [p2| , can be written, up to terms of order a 2 s , as: 



r^up (x) 



Gla 



4~ m b,pole m b -MS 



32tt 1 

QQ*b) 



(m b ) \V tb V t *\ 



* I 2 



9(1- x) - ^A x a s (n b ) log 2 (1 - x) + B x a s (ji b ) log (1 - x) 



+ p(x) 



(124) 



■Oi 



One then takes a derivative with respect to x. The requirement that the spectrum is non-singular at x = 1 
or equivalently, that the total rate is correctly reproduced^ — transforms the distributions in plus distributions 
and one obtains eq. (116). 



4 Conclusions 



We have performed factorization and threshold resummation of the most general distribution in the semileptonic 
b — > u decay. This has been achieved with a proper choice of the kinematic variables, 



2E X 

m B 4E X 



and £= (125) 



which disentangle the two different logarithmic structures occurring in the process: 

1): \ogw and 2): log(l-z). (126) 

The first structure has been related to the infinite mass limit of the b quark, 

tub — * oo, (127) 

while the second one has been related to the infinite energy limit of the final hadronic system, 

E x -> oo. (128) 

The long-distance effects — both perturbative and non-perturbative — have been relegated to a universal factor 
/ (z) , depending only on z. This function is related to the shape function ip (k + ) by a short-distance coefficient 
function. In essence, in our approach, the hard scale of the process is the hadronic energy Ex in the B rest 
frame, and not the B meson mass tub- 

We have presented simple analytical expressions for a few distributions, resummed to leading logarithmic 
accuracy and we have shown that the shape function can be directly determined by measuring the distribution 
in z or in z and w. 



6 We do not consider here the problem of the infrared singularities for x — > 0, related to soft-photon emission. 
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We have computed the resummed hadron energy spectrum, which exhibits a "Sudakov shoulder" in the 
point 

Ex = ^ (129) 

i.e. in the middle of the allowed kinematical domain (0 < Ex < to_b) . The spectrum very close to this point is 
non-perturbative and is proportional to the shape function. This implies that, at least in principle, an accurate 
measure in this region may lead to another independent determination of the shape function. It is a non-trivial 
fact that the same long-distance effects appear at the boundary of the phase space in the z-distribution while 
they appear inside the allowed kinematical domain in the energy spectrum. The computation of other resummed 
distributions in LO or NLO by integrating eq. ( |90| ) is straightforward. A crucial point in our analysis is that 
there is a single source of large logarithms in any semi-inclusive distribution in heavy-flavour decay. 



Finally, we have presented an improved formula for the photon spectrum in the rare decay (25), which takes 
into account soft-gluon resummation and the effects of the subleading operators. It has been explicitly shown 
that the same function / (z) factorizes the long-distance effects in the semileptonic and in the rare decay. 

We believe that our formalism sets a rather general and rigorous scheme for the separation of perturbative 
from non-perturbative effects in semi-inclusive heavy-flavour decays and allows for a simple analysis of the 
experimental data of many different distributions. 
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